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SYMBOLS 


To create a degree of order, the following rules have been adopted in the assignment of mathematical 
symbols to the various quantities. 

• All variables are designated by lower case symbols (Roman or Greek). 

• With the exception of the acceleration due to gravity, g, all constants are designated by upper 
case symbols (Roman or Greek). 

• A hat over a quantity (e.g., x) denotes that the quantity has been normalized. Quantities 
associated with the acquiring curve are normalized by the acquiring curve length, L, and quantities 
associated with the blending curve are normalized by the true circle radius, &i 

• With the exception of (adopted to conform with ref. 2), a primed quantity (e.g., x ') denotes 
one that has been differentiated with respect to the range, d. 

• A dot over a quantity (e.g., x) denotes differentiation with respect to time. 

• Vector quantities are shown boldface (e.g., r). 

The symbols are divided into the following four categories: 

1 Input variables (either from sensors or precomputed and stored). 

2. Constant inputs (set by the pilot or frozen at flightpath select). 

3. Derived variables used in the description of the flightpath and guidance laws. 

4. Derived constants used in the description of the flightpath and guidance laws. 

Input Variables 


h 

aircraft altitude, feet (ft) 


h 

vertical velocity of aircraft’s center of gravity (e.g.), feet per second (ft/sec) 

vt 

aircraft groundspeed, ft/sec 


vt 

aircraft acceleration along ground track, ft/sec 2 


Vxi Vy 

horizontal longitudinal and lateral speeds, ft/sec 


V x ,Vy 

horizontal longitudinal and lateral accelerations, ft/sec 2 

X 

normalized blending-curve ordinate 


%a <> Ua 

aircraft coordinates in the landing-pad coordinate 

system, ft 

St 

throttle position, degrees (deg) 


Stc 

throttle compensation parameter, deg 


Sxi Sy 

longitudinal and lateral pilot control inputs, deg 


C 

direction of ghost aircraft relative to real aircraft, 

deg 

6 

aircraft pitch angle, deg 



SB -W 4$^- 
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<7 W aircraft vertical damping parameter, 1/sec 

a x ,cr y aircraft horizontal translational damping parameters, 1/sec 

4> aircraft roll angle, deg 

Input Constants 

Cq, pitch and roll rate sensitivity to pilot control inputs 

D c nominal length of final flightpath segment, ft 

Df threshold range, ft 

D x , D y aircraft pitch and roll damping constants, 1/sec 

Eq heading of head-up-display (HUD) flightpath symbol at flightpath select, deg 

g acceleration due to gravity, ft/sec 2 

H hover altitude, ft 

H r i aircraft altitude at guidance select, ft 

K a b ratio of lengths of acquiring curve and initial straight segment 

K e flightpath-symbol lateral-drive scale factor 

Kg x , Ks y control-input gains, feet per inch second (ft/in. sec) 

L\ normalized distance along acquiring curve, measured along x-axis, where the angle of 

bank changes sign (fig. 6) 

Mf,Hf azimuth and elevation of a fixed point in space far ahead of the aircraft, deg 

Rl(min) nominal minimum radius of horizontal flightpath circles, ft 

R radius, ft 

Rv radius of vertical flightpath circle, ft 

S x , S y touchdown or station-keeping point coordinates in the landing-pad reference frame, ft 
2*3 control-input washout time constant, sec 

Ty horizontal-tracking time constant, sec 

V x f desired value of v x at threshold range, Df, ft/sec 

X u effective value of sign reversed slope of speed-versus-range line of reconfigured aircraft, 

1/sec 

r approach flightpath angle, deg 

A D ac minimum desirable distance along curved segment of flightpath, ft 

A T g general ghost lead time, sec 

Nq (mi n ) minimum permitted value of % deg 

Tp normalized maximum length of portion of ordinate of blending curve replacing a portion 

of a straight segment 

T normalized length of portion of ordinate of blending curve replacing a portion of a 

straight segment 

W landing-pad heading, deg 

^ f Heading of final flightpath segment, deg 

£l 0 translational velocity response frequency, radians/sec (rad/sec) 
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Derived Variables 


bx ) by 

d 

dy* 

lc 

he 

kg 

hf 

i 

k W g 

k 8 t 

lac 

r a 

r 9 

n 

m 

r llxi r lly 
S 

s{x) 
s(x p ) 


Vxi'by 
Vxci Vyc 
Vxr 

x,y 

x f>yf 

x,y 
x a ,Va 
Xp, y p 


x\ai yia 
X2a , V2 a 


It 

6h 

SXg 

fiXgh 

6y 

Sy 


coordinates of point B in landing-pad system prior to flightpath select, ft 

range to touchdown or station-keeping point measured along synthesized flightpath, ft 

range at which aircraft should be reconfigured, ft 

distance from aircraft to center of flightpath circle, ft 

quickened vertical velocity, ft/sec 

altitude of ghost aircraft, ft 

altitude of reference flightpath at aircraft’s range, ft 

indicates right (i = +1) or left (i = — 1) flightpath circle prior to flightpath select 
ghost lead blending gain 
throttle washout gain, ft/sec deg 

distance of aircraft from tangent point C (fig. A-3) of circle with final straight segment, ft 
vector position of real aircraft, ft 
vector position of ghost aircraft, ft 
instantaneous circle radius, ft 

radius of appropriate circle prior to flightpath select, ft 

coordinates of appropriate circle in landing-pad reference frame prior to flightpath select, 
ft 

Laplace transform variable, 1/sec 

nondimensional distance along blending curve from origin (fig. 5) 
nondimensional distance along blending curve from origin to foot of perpendicular from 
aircraft to blending curve (fig. 5) 

smoothed estimates of horizontal translational accelerations, ft/sec 2 
commanded horizontal translational velocities, ft/sec 
reference longitudinal acceleration, ft/sec 2 
coordinates of acquiring curve, ft 

HUD coordinates of a fixed point in space with azimuth and elevation of Mf and Hf, 
respectively, deg 

nondimensional coordinates of blending curve (fig. 5) or acquiring curve (fig. 6) 
nondimensional coordinates of aircraft in x, y axes system (fig. 5) 
nondimensional coordinates of foot of perpendicular from aircraft to blending curve in 
x, y axes system (fig. 5) 

coordinates of aircraft in an axis system with origin at center of flightpath circle and 
x-axis parallel to final flightpath segment, ft 
coordinates of aircraft in an axis system with origin at point B (fig. B-l) and x-axis 
parallel to initial flightpath segment, or with origin at point S and x-axis along final 
flightpath segment, ft 

reference flightpath angle at aircraft’s position, deg 

altitude error relative to reference glide slope, ft 

general ghost lead distance, ft 

ghost lead distance for vertical tracking, ft 

lateral distance of aircraft from reference flightpath track, ft 

nondimensional perpendicular distance from aircraft to blending curve (fig. 5) 
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Sipi 

£ q 

i 

Hg,r)g 

V 

^0 

?° 

% 

p 

V 

4 } g 

4 ) max 

4>ac 


A(i),A( 2 ) 

B (1)’ B (2) 
Bx -> By 

C (1)' C (2) 

C X ,Cy 


Dai 

D bi 

D b 2 

D c 

Dei 

Dc 2 

B j 

B h 

B J2 

E 

Eg 

F 

I 

L 


L i 


heading of initial flightpath segment relative to x-axis of landing pad prior to flightpath 
select, deg 

heading of flightpath relative to heading of aircraft prior to flightpath select, deg 
nondimensional incremental distance between x p and x a 
azimuth and elevation of ghost aircraft, deg 

angle subtended at center of reference flightpath circle by remaining arc of circular 
flightpath segment, deg 

nominal angle subtended at center of reference flightpath circle by entire circular 
flightpath segment, deg 

angle between line AC (fig. A-2) and final straight segment, deg 
limited value of Pq, deg 

radius of curvature of reference flightpath track, ft 

nondimensional parameter defining length of blending curve prior to flightpath select 
ghost-aircraft roll angle, deg 

minimum value of maximum angle of bank required during flight along acquiring curve, 
deg 

heading of line AC (fig. A-2), deg 

heading of tangent to reference flightpath track at aircraft’s position, deg 

Derived Constants 

constants 

constants 

coordinates of point B in landing-pad system after flightpath select, ft 
constants 

coordinates of point C in landing-pad system, ft 

range at point A, ft 

range at point A\ (fig. B-l), ft 

range at point B\ (fig. B-l), ft 

range at point B 2 (fig. B-l), ft 

range at point C (fig. B-l), ft 

range at point C\ (fig. B-l), ft 

range at point C 2 (fig. B-l), ft 

range at point J (fig. B-2), ft 

range at point J\ (fig. B-2), ft 

range at point J 2 (fig. B-2), ft 

angle subtended by blending curve at center of reference flightpath circle, deg 
value of e q 4- after flightpath select, deg 

nondimensional quantity used in evaluation of X 

indicates right ( I — +1) or left (/ = —1) flightpath circle after flightpath select 
length of acquiring curve measured along x-axis (fig. 6), ft 

distance along acquiring curve, measured along x-axis, where angle of bank changes sign 
(fig. 6), ft 
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R ll 

R llx^ R lly 

R b R ll 

R l(max) 

R l 

p(<) 

jrt Z(mm)x 

»(*) 

' Lh l(min)y 

V' 

v w 

A 

A Vi 
A 

H r i 

N 0 

N 0 


X 


^ -vi 


HUD 

IMC 

VMC 

VMS 

VSRA 


distances from aircraft to centers of flightpath circles, ft 
nominal radius of flightpath circle after flightpath select, ft 
coordinates of center of flightpath circle, ft 
true radius of reference flightpath circle, ft 
nominal maximum radius of horizontal flightpath circles, ft 

nondimensional nominal radius of reference flightpath circle 
x-coordinates of minimum-radius circles in landing-pad reference frame, ft 

y-coordinates of minimum-radius circles in landing-pad reference frame, ft 
effective wind speed, ft/sec 

heading of final flightpath segment relative to x-axis of landing pad, deg 
angle between line from aircraft to center of reference flightpath circle and x-axis of 
landing pad, deg 

heading of initial flightpath segment relative to x-axis of landing pad after flightpath 
select, deg 

aircraft altitude at flightpath select, ft 

angle subtended at center of nominal reference flightpath circle by complete arc of 
flightpath’s circular segment, deg 

angle subtended at center of true reference flightpath circle by complete arc of 
flightpath’s circular segment, deg 
nondimensional coefficient of blending curve 
nondimensional ordinate length of entire blending curve 

nondimensional ordinate length of portion of blending curve lying within arc of nominal 
flightpath reference circle 
heading of aircraft at flightpath select, deg 
heading of initial, straight flightpath segment, deg 

heading of velocity vector relative to initial, straight flightpath segment at flightpath 
select, deg 

Acronyms 

head-up display 

instrument meteorological conditions 
visual meteorological conditions 
vertical motion simulator 
V/STOL systems research aircraft 


i 


Vlll 



SUMMARY 


A piloted simulation test was conducted on the Ames Research Center’s vertical motion simulator 
(VMS) in support of the Phase HA flight test of NASA’s V/STOL systems research aircraft (VSRA). 
During the simulation several problems were found with the head-up display (HUD) symbol drive laws 
and the flightpath synthesis. These problems and the solutions devised to solve them are described in 
this report. Most of the resulting HUD drive-law changes were implemented during the simulation and 
their effectiveness was verified. Subsequently both the HUD symbol drive-law and flightpath-synthesis 
changes were implemented in the VSRA and tested successfully in the Phase HA flight tests. 

INTRODUCTION 

The Phase HA flight tests of NASA’s VSRA, a YAV-8B Harrier, were performed in late 1991. 
For the Phase HA tests, the aircraft was equipped with a programmable HUD, flight computers, and 
sensors to provide position, velocity, and acceleration data. The aircraft control system used in the tests 
was a standard, inservice Harrier. The purpose of the flight tests was to evaluate an advanced HUD 
format, HUD symbol drive laws, and a flightpath-synthesis technique aimed at reducing the visibility 
minima for landing. The HUD format, shown in figure 1, is described in detail in reference 1. The 
flightpath-synthesis technique and the HUD symbol drive laws, with the exception of an exponential 
speed-guidance law, are also described in reference 1. The exponential speed-guidance law is one that is 
suitable for the type of control system used in the tests; it is described in reference 2. Because the aircraft 
has a single seat, safety considerations dictated that the flight tests be performed in visual meteorological 
conditions (VMC). In an attempt to extend the flight test results to instrument meteorological conditions 
(IMC), a piloted, moving base simulation was conducted on the Ames Research Center’s VMS; the 
visibility conditions were varied from VMC to “zero-zero.” 

The simulation was conducted just prior to the Phase DA VSRA flight tests. During the simulation, 
a variety of problems with the HUD drive laws and flightpath synthesis were uncovered. The purpose 
of this report is to describe these problems and the changes devised to solve them. Those changes 
whose implementation required relatively minor code modifications were made and tested during the 
simulation. All the changes, including those whose implementation required more extensive code 
modifications, were tested during the flight tests. 

The report addresses seven distinct changes to the material contained in references 1 and 2, as 
follows: 

1 . A fundamental error in reference 1 is identified and corrected. This error concerns the way in 
which the flightpath and ghost symbols (fig. 1) were scaled laterally. 

2. Pilots noted that, when following the vertical guidance during the deceleration phase after 
reconfiguration (nozzle drop), the flightpath was invariably oelow the glide-slope reference (fig. 1). The 
reason for this vertical error and a correction for it are identified. 

3. The speed-guidance law given in reference 2 did not have sufficient flexibility to modify ad- 
equately the speed characteristics in the final few hundred feet of range prior to reaching the hover 
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Approach HUD 
Distance to Landing 



Hover HUD 

Figure 1 HUD format. 


point. This fact became clear when an experienced Harrier pilot requested that the final approach be 
quickened. The problem is shown to stem from insufficient consideration of the threshold speed and 
range in the speed-guidance law. An improved speed-guidance law is derived. 

4. Recent research (ref. 3) has provided a more rational approach to deriving laws for the velocity- 
predictor ball (fig. 1) than that used to produce the law given in reference 1. Apart from improved 
dynamics, this new approach permits a simple means of modifying the ball response independently of 
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the damping factor and is easier to adjust to pilot requirements. A modified version of the law given in 
reference 3 is derived to replace that given in reference 1 

5. Even when the ghost aircraft was tracked precisely, significant lateral flightpath errors were 
noted at turn entrances and exits, most notably in the latter. It is shown that this anomaly is caused by 
discontinuities in the lateral curvature of the reference flightpath at points where straight and circular 
segments meet. The proposed solution given herein is to remove these discontinuities by providing a 
“blending curve” in the region where the straight and circular flightpath segments join. This blending 
curve is designed to provide a reference flightpath with a continuous curvature. 

6. It has been noted that if, at the instant of flightpath select, the heading of the aircraft’s track 
differs significantly from that of the reference flightpath, then the required turn to acquire the reference 
flightpath can be large enough to cause the ghost-aircraft symbol on the HUD to saturate. If, under these 
circumstances, the pilot is using the HUD for guidance, the ghost gives the direction of the turn but 
not its magnitude. This lack of information can cause the pilot either to overcontrol or to undercontrol. 
Overcoming this problem involves defining an “acquiring curve” that provides explicit guidance for the 
“S” turn to the first straight flightpath segment. 

7. In reference 1, no explicit vertical guidance was provided other than a glidepath. The pilot was 
expected to fly the aircraft until glidepath intercept occurred, as indicated by the descent of the ghost 
symbol on the HUD. This indication was too abrupt and it resulted in significant glide-slope overshoot. 
To remedy this problem, the reference flightpath now maintains a constant altitude equal to that of the 
aircraft at the instant of flightpath select and is gradually blended into the glide slope. This modification 
is included in the appendixes. 

Comprehensive sets of equations for the flightpath synthesis and guidance incorporating the adopted 
solutions to the problems are given in appendixes A and B. Results of the simulation and flight tests 
are given in reference 4. 


HUD LATERAL SCALING 

During the simulation, it was noted by one of the pilots that when performing steady, coordinated 
turns in calm conditions, both the ghost and flightpath symbols were displaced laterally from the HUD 
centerline. Although the flightpath tracking task appeared to be satisfactory, it was clear that this lateral 
displacement should not occur, since it erroneously signified a sideslip. The problem resulted from the 
manner in which both the flightpath and the ghost symbols were scaled laterally. Lateral scaling is 
needed because of HUD field-of-view limitations. The HUD projection equipment in the VSRA has a 
usable lateral field-of-view of about ±1 deg. Since, at very low speed, in a crosswind, the aircraft yaw 
relative to the flight direction can easily exceed 20 deg, it is clear that some degree of lateral scaling is 
needed. In both the simulation and the flight test, the lateral motion of the flightpath and the ghost was 
reduced to 30% of the full value. In the notation of reference 1, this scaling corresponds to a value of 
K e of 0.3. In reference 1, this scaling factor was applied to the lateral motion of both the ghost and the 
flightpath symbols when defined in Earth-fixed axes prior to the transformation to HUD axes. Clearly, 
this procedure is in error since the intention was to scale the lateral motion in HUD axes. It follows that 
the correct approach is to reverse the procedure by applying the transformation first and then applying 


j 


3 



the lateral scaling factor. The appropriate corrections to the symbol drive laws of reference 1 are to 
delete the K e factors from equations (15) and (23), which define the lateral positions of the flightpath 
and the ghost symbols in Earth-fixed axes, and to multiply by K e the lateral positions in HUD axes 
given by equations (17) and (25). These corrections were made for both the VMS simulation and the 
Phase HA VSRA flight tests with satisfactory results, at least for the approach task. 

Although the new scaling procedure just described is fundamentally sound, further experiments 
have shown that it results in another problem, noted particularly during large changes of roll angle. In 
order to see the nature of this problem, consider the path traced out on the HUD by a point fixed in 
space and far ahead of the aircraft, when the aircraft is rolled at constant pitch angle through 360 deg. 
Let the azimuth and elevation of the fixed point be Mf and Hf, respectively. The unsealed HUD 
coordinates of this point, Xf and yf, are given by 

Xf = Mjcoscf) - (Hf - 6) sin0 
yf = Mf sin 04- (Hf - 0)cos0 
where 6 and 0 are the aircraft pitch and roll angles, respectively. 

Eliminating 0 between these equations gives 

x 2 f + y 2 f = M 2 f + (H f - 9) 2 

from which it follows that, if Mf, Hf, and 6 are all constant, the path traced out on the HUD by the 
fixed point is a circle of radius \jMj + (H'j — 0 2 ). Suppose now that the lateral position Mf is scaled 
by a factor K e , less than unity, as was done in refe rence 1 The above e quations show that the path 

of the point is still circular, but of reduced radius ^K e 2 Mf + (Hf — 9 2 ) Alternatively, if the lateral 

HUD coordinate Xf, rather than Mf, is scaled by K e , then the path of the point on the HUD is the 
ellipse 

x 2 f K 2 + y 2 f = M 2 f + (Hf-e ) 2 

oriented so that its major axis is vertical and the ratio of its minor to major axes is K e . 

It follows from the above analysis that, although the method of lateral scaling given in reference 1 
produces incorrect angles of sideslip, the path of a fixed point, which could be either the ghost symbol or 
the flightpath symbol, is at least geometrically similar to that which would be observed if the HUD were 
unsealed. On the other hand, the new method of scaling, although not showing false sideslip angles, 
produces a very distorted representation of the path of a fixed point. It turns out that this distortion, 
although not too noticeable in a precise tracking task where only moderate rates of roll are required, 
does cause problems in circumstances where large rates of roll and large roll angles are required. The 
problem occurs because the fixed point appears to move much more vertically than laterally, giving the 
pilot a strong visual cue that a vertical-control correction is required. 

Clearly there is no way to disguise the fact that the lateral motion is being scaled. The scaling will 
produce an anomaly of one form or another. The best that can be done is to perform the scaling in a 
way that is the least dangerous. A possible solution is to remove all scaling at speeds high enough that 
crab angles will not normally exceed the lateral field-of-view of the HUD. The new form of scaling can 
be introduced at the instant when the terminal guidance is selected, since it is only beyond this point 
that crab angles can become large. 
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ALTITUDE ERROR DURING APPROACH DECELERATION 


As explained in detail in reference 1, the HUD concept used in the VMS simulation and Phase HA 
VSRA flight tests employs a ghost aircraft to represent a ficticious aircraft flying ahead of the real aircraft 
and performing the approach task perfectly. To replicate this desired approach, the pilot maneuvers the 
real aircraft so that the flightpath symbol (velocity vector) on the HUD is coincident with the ghost 
aircraft. 


During the early phase of the VMS simulation, it was found that during the approach deceleration 
after reconfiguration (nozzle drop) the aircraft showed a pronounced tendency to fall below the desired 
flightpath. The maximum altitude error was typically about 30 ft. This altitude error occurred even 
though the pilot maintained the flightpath symbol precisely on the ghost symbol. The altitude error was 
qualitatively evident to the pilot since, even though the two symbols were coincident, both were above 
the glide-slope reference line on the HUD (fig. 1). Surprisingly, this significant inconsistency on the 
HUD passed unnoticed in the simulation tests of reference 2. 


This inconsistency problem was caused by the vertical quickening introduced into the flightpath- 
symbol vertical-drive law. This quickening was produced by augmenting the vertical velocity with 
“washed-out” throttle position. The effective vertical velocity used to calculate the vertical position of 
the flightpath symbol on the HUD (vertical flightpath angle) is given in reference 1 by 


he — h + 


S + &W 


( 1 ) 


where 

h 

vertical velocity 

he 

quickened vertical velocity 

St 

throttle position 

k 6 t 

throttle washout gain 

s 

Laplace transform variable 

G W 

aircraft vertical damping parameter 


The reason for using quickened vertical velocity in the flightpath-symbol drive law is to make easier 
the pilot’s task of maintaining the desired flightpath, and it undoubtedly achieves this goal. However, 

h c = h in the long term only if the average deviation of St from some nominal constant is zero in the 
long term. This condition is met in conventional, constant-speed approaches. In decelerating approaches, 
the throttle must be advanced continuously to compensate for the gradual loss of aerodynamic lift as 
the airspeed decreases. In this case, the rate of change of St (i.e., sSt in eq (1)) is positive for as long 

as 30 sec, during which time h c ^ h. This “standoff” error is exaggerated by the long washout time 

constant of about 10 sec. Since the HUD flightpath symbol depends on h c rather than h, it 

does not provide the pilot with the correct flightpath angle. More specifically, during the deceleration, 

h c > h and the HUD-indicated flightpath angle (measured positive when climbing) is greater than the 
true flightpath angle by as much as 2 deg. It follows that, if the pilot maintains the flightpath symbol on 
the ghost, then the true flightpath angle will be less (more negative) than the indicated flightpath angle 
and the aircraft will settle below the desired flightpath. 
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To overcome this difficulty, equation (1) was modified to be 


he — h + 


htfjSt - he) 

S + (Tyj 


( 2 ) 


where 6tc is the throttle angle required to maintain the desired flightpath angle during the deceleration. In 
an effort to determine a suitable form for 6t c , a pilot was requested to perform a series of decelerating 
approaches, in various winds, using the original flightpath-symbol vertical-drive law. The pilot was 
instructed to use the HUD to keep the ghost on the glide-slope reference line, even though this meant 
keeping the flightpath and ghost symbols noncoincident. The resulting variation of throttle angle with 
airspeed are shown in figure 2. It can be seen that the correlation of the data for the various approaches 
is sufficiently good that the data can be approximated by the single curve of figure 3, This curve was 
used to define the value of 6f C in equation (2). The modification shown in equation (2) overcame the 
problem and turned out to be satisfactory over a range of flightpath angles from -2 deg to -5 deg 
(ref. 4). 



Figure 2. Variation of throttle angle with filtered airspeed during deceleration. 





Figure 3. Throttle compensation used in flightpath-symbol vertical-drive law 


SPEED-GUIDANCE LAW MODIFICATION 


The speed-guidance law given in reference 2 consists of two parts, The first part specifies the 
reconfiguration range, d r , where the engine nozzle control is moved to the hover position. It is defined 
to be the range, d, at which the inequality 


is first satisfied, where 



(3) 


d aircraft range from the selected station-keeping point 

d r reconfiguration range 

Df threshold range 

v x horizontal longitudinal speed relative to the station-keeping point 
V x f desired value of v x at the threshold range, Df 

Vlj effective wind speed (ref. 2) 

X u effective value of the sign reversed slope of the speed-versus-range line of the reconfigured 
aircraft (ref. 2) 


The second part of the speed-guidance law specifies that the pilot produce a deceleration, after 
reconfiguration, given by 


Vxr — — 


v: 


x 


d 


(4) 
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where v xr is the reference horizontal longitudinal acceleration. 


In the piloted simulation Df and V x f were initially set at 200 ft and 30 ft/sec, respectively. 
However, an experienced Harrier pilot, flying in VMC, judged the final approach from the threshold to 
be too slow. Attempts were made to speed up the final approach to the hover point by changing Df and 
V x f, but these changes were unsatisfactory. Changing Df and V x f changed only the reconfiguration 
range (eq. (3)), and the guidance law (eq. (4)) always tended to cancel the effect of this change so that 
the change in the speed of the final approach was only minor. The problem with the overall guidance 
scheme given above is that the deceleration law is incompatible with the reconfiguration range. The 
exact deceleration law is given by equation (5) of reference 2, namely 

tto- = *«(% + 0 (5) 

where the horizontal longitudinal acceleration v x has been treated as the reference value v xr . 


The problem with this law is that it is strongly dependent on X u , which is dependent on the 
average angle of attack needed during the approach, which, in turn, is dependent on the reconfiguration 
range. Therefore, unless the pilot moves the engine nozzle to the hover position at precisely the 
reconfiguration range, the guidance law will not provide guidance to the desired threshold conditions. 
In essence, equation (5) is not a feedback law and it does not have any continuous acceleration-correction 
properties. What is required is a version of equation (5) that specifically contains the range d and the 
threshold parameters Df and V x f. This version can be achieved by first noting that, during a perfect 
approach in which the angle of attack is constant and is the correct value, the relationship between 
speed and range is given by 


d=- 


f). 

and this equation may be rearranged to give 


/X u + L 


( 6 ) 


X u = 



which, when substituted into equation (5), gives 


(7) 


v x - 



( 8 ) 


Equation (8) is a more exact version of equation (4); it shows clearly that implicit in equation (4) is the 
assumption that = D f = V x f = 0. However, the results of both the simulation on the VMS and the 
simulation reported in reference 2 indicate that, when equation (4) is used for speed guidance, nonzero 
values of were compensated for quite well. Therefore, it is reasonable to presume that the same 
situation would be true if equation (8) were used for speed guidance. Setting = 0 in equation (8) 
gives 


Vxr — 


d- 




(9) 


X 


8 



This equation is the guidance law that emerged from the VMS simulation as being satisfactory (ref. 4). 
Although it depends on X u , the dependency is not so strong that the speed guidance is significantly 
affected by errors in the reconfiguration range. 


VELOCITY-PREDICTOR-BALL LAW 


The velocity-predictor-ball law used in the VMS simulation and Phase HA VSRA flight tests 
differed somewhat from that given in reference 1. In reference 1 the commanded horizontal velocities 
v X c and v yc are given by the expressions 


- .. , m 7 . K 6 x T 2 sS x 

V X C V X + TyV X 'J'nS + 1 

(10) 

K$ T2s6 y 

Vyc = Vy+ TyVy + T ^ g+1 

(11) 


where 

v X c, Vyc commanded horizontal translational velocities 

v x , v y horizontal translational speeds 

ti x , v y smoothed estimates of the horizontal translational accelerations 

8 X , 8 y longitudinal and lateral pilot control inputs 

K § x , Kg y control-input gains 

T 2 control-input washout time constant 

Ty horizontal-tracking time constant 

This law was deduced in an ad hoc manner with the gains Kg x ,Kg ,T 2 , and Ty set, experimentally, in 
piloted simulations. 


A more rationally designed predictor-ball law was derived in reference 3 and tested on a variable- 
stability helicopter. One of the attractive features of this law is that its principal dynamic character- 
istic (break frequency) can be varied through a single frequency parameter, Oo, while maintaining a 
closed-loop, third-order, binomial transfer function for the horizontal velocity response to a commanded 
predictor-ball position. Since all the gains are functions of Qq, it is easy to change the speed of response 
of the aircraft to a given commanded ball position. A variation of the approach used in reference 3 
was used in the VMS simulation and Phase HA VSRA flight tests. The third-order, binomial transfer 
function was chosen to represent the response of the aircraft horizontal position rather than the hori- 
zontal velocity. The steps in the derivation of the corresponding predictor-ball law closely follow those 
outlined in reference 3. The resulting law is 


V XC 


Vi 


yc 


2 ^ gD x CgS(s + 40o)^ar 

flo X 4f2o 3 (s + 0*){« + D x ) 

(12) 

2 - 9DyC^s(s -1- 4f2o)^y 

Hq ^ + &y){s + Dy) 

(13) 
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where 

O 0 translational velocity response frequency 

g acceleration due to gravity 

D x , D y aircraft pitch and roll damping constants 

Cq, pitch and roll rate sensitivity to pilot control inputs 

a x , <j y aircraft horizontal translational damping constants 

AIRCRAFT/GHOST FLIGHTPATH GEOMETRY 

The kinematical implications of a flightpath guidance technique that uses a ghost aircraft flying 
ahead of the real aircraft should be examined before we consider a method of reducing flightpath errors 
at turn entrances and exits (problem number 5 in the Introduction). Consider the situation depicted in 
figure 4. The real aircraft is being maintained on a predetermined flightpath (solid line) such that its 
velocity vector always points to a ghost aircraft located on an adjacent flightpath (dashed line). Let the 
flightpath of the aircraft be given by 

r a = r aid) (14) 

where r a is the vector position of the real aircraft measured from an origin 0 and d is the aircraft’s 
range measured from a fixed point A on the flightpath (fig. 4). The unit tangent vector to the aircraft’s 
flightpath at the position of the aircraft is r^, where the prime denotes differentiation with respect to the 
range, d. It follows from the geometry of figure 4 that the position of the ghost aircraft, r g , is given by 



Figure 4. Aircraft/ghost flightpath geometry. 
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where 8x g is the distance from the ghost aircraft to the real aircraft. The unit tangent vector to the 
ghost aircraft’s flightpath, r' g , is obtained by differentiating equation (15) with respect to d; thus, 

4 = 4 + 8x g r" + 8x g r' a (16) 

The direction of the ghost aircraft relative to the real aircraft, ( (fig. 4), is given by the inner product 
of r'g and r' a , as follows. 


4'4 = 1*5 1 1 ""a I cos C - (r a + 6x 9 r a + &44)'4 


(17) 


r ' .j/ — ij*^ p == x 
l a *a l^ai x 

in 


Since is a unit vector, 

K a*a I ‘a 

Differentiating equation (18) with respect to d results 

r' *r" = 0 
l a l a u 

Applying equations (18) and (19) to equation (17) shows that 

Ir^lcosC = 1 + 8x' g 

But it follows from equation (16), with equations (18) and (19), that 

\*'g\ 2 = [4(1 + 8x 'g) + %4'l [4(1 + Sx 'g) + Sx g*a\ 

= 44(1 + Sx’gf + 8x g 2 4-4 + 28x g (l + < 54 ) 4=4 

= (i + K) 2 + (t 2 ) 2 

where p is the radius of curvature of the real aircraft’s flightpath at the position of the aircraft. 
Specifically, 

def 1 

P = t 


It follows from equations (20) and (21) that 

cos C, = 




1 + 


6Xr 


(18) 

(19) 

( 20 ) 

(21) 

( 22 ) 

(23) 


The term 8x g expressed as a function of the aircraft’s kinematical variables is: 

. j d8x g d8x g dt 8x g 

dd dt dd vt 


(24) 


where vt is the groundspeed of the real aircraft. If, as is usual, the ghost lead distance 8x g is made 
proportional to vt, that is, 

8Xg = VfATg (25) 


j 
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where A T g is the (constant) ghost lead time, then differentiating equation (25) with respect to time and 
substituting for 6x g in equation (24) gives the following final expression for 6x' g : 


y vt 

where vt is the aircraft’s acceleration along the ground track. 


(26) 


It is instructive to evaluate the implications of equations (23) and (26) by means of two examples. 


Example 1 

An aircraft, traveling at a constant groundspeed, vt, of 200 ft/sec (vt = 0) on a circular flightpath of radius, 
p, of 4,000 ft, is just about to turn onto a straight flightpath segment (p = oo). The lead time of the ghost 
aircraft, A T g , is 10 sec. It follows from equation (25) that the lead distance of the ghost aircraft, 8x g , is 2,000 ft 
and from equation (24) that 6x' g = 0. Just prior to and just after exiting the turn, ( = 26.6 deg and C = 0 deg, 
respectively (eq. (23)). Therefore, in exiting the turn, the ghost-aircraft flightpath instantly changes direction by 
26.6 deg relative to the real-aircraft flightpath. The pilot perceives a sudden change of direction of this magnitude 
as an undesirable heading transient despite the measure taken in reference 1 to smooth it out by introducing a 
turn-exit blend distance. 


Example 2. 

As in example 1, the aircraft is traveling at a constant groundspeed, vt, of 200 ft/sec (vt = 0) on a circular 
flightpath of radius, p, of 4,000 ft and a ghost-aircraft lead time, AT g , of 10 sec. However, in this example, the 
aircraft is suddenly decelerated at 8 ft/sec (i>t = —8 ft/sec) while still turning. A deceleration of this magnitude 
is produced on a Harrier, typically, when the engine nozzle is moved to the hover position. At the instant of 
deceleration, the value of 8x g changes from zero to —0.4 (eq. 24)) and the measure of ( changes from 26.6 deg 
to 39.8 deg (eq. (23)). This 13.2-deg change of £ is probably not discernible to the pilot. However, if future 
flightpath requirements are for smaller radii circles to be flown very accurately, then the effect could become 
more important. If, for example, the radius were to be 2,000 ft instead of 4,000 ft, then the change of C would 
be a much more significant 32.4 deg. 


BLENDING-CURVE GEOMETRY 


The previous analysis of the aircraft/ghost flightpath geometry indicates that it is desirable to have 
a synthesized flightpath with continuous curvature. One way to achieve this continuity, while retaining 
much of the synthesis technique of reference 1, is developed in this section. 

The initial process of synthesizing a flightpath trajectory involves connecting initial and final straight 
segments by a circular segment with a given nominal radius, This procedure produces a flightpath 
whose curvature is continuous except at the two points of tangency of the straight and circular segments. 
The point of tangency of one of the straight segments and the circle is shown in figure 5 at point 2. The 
approach adopted here to remove the curvature discontinuity is to modify the geometry of the flightpath 
locally by means of a “blending curve.” The process of adding a blending curve between the straight and 
circular segments starts by defining a second circular segment with the same center as the original but 
of slightly smaller radius, R[, and connecting this new circular segment and the straight segment with 
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Figure 5. Flightpath transition-curve geometry. 

a cubic blending curve, 01 (fig. 5). This blending curve is defined to have, at point 0, the same slope 
and zero curvature as the straight segment and, at point 1, the same slope and curvature as the circle. 
Since the blending curve is a cubic, its curvature is approximately linear along its length, varying from 
zero at point 0 to l/Ri at point 1. The complete synthesized flightpath consists of a straight segment, 
a blending curve, a circle of radius R\, a second blending curve, and a final straight segment. All the 
dimensions in figure 5 have been normalized by the radius R\ so that the true flightpath circle has unit 
radius. Normalization is uniformly indicated by a hat over the symbol. 

The blending curve has the form 

y = Ex 3 (27) 

As shown in figure 5, the origin of the blending curve is defined to be at a normalized distance T from the 
intersection of the straight segment and the nominal flightpath circle (point 2 of fig. 5). Equation (27) 
satisfies the required conditions of continuity, slope, and curvature at point 0. The task now is to 
determine the value of 5 consistent with the required continuity, slope, and curvature at point 1. 
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The equation of the true flightpath circle in the coordinate system shown in figure 5 is 


(:r-f f + {y-Rj ) 2 (28) 

The following relationships express the required continuity of position, slope, and curvature (second 
derivative) at point 1* 


position (X - T) 2 + (EX 3 -Rif = 1 


(29) 


PX 2 = - 


slope 3EX 


curvature . 6EX = 


From equations (30) and (31), 


1 + 


/ X-f \ 
\EX 3 - R\) 

EX 3 - Ri 


( x-r \ 2 

l ex^M 


(30) 


(31) 


(32) 

X X-T 

Substituting for (EX 3 — Rif from equation (29) into equation (30) and simplifying gives the following 
cubic equation: 


where the quantity X is defined as 


^3 X T „ 
X - — + — = 0 
2 2 


X = f X - f 


The appropriate real, positive solution of equation (33) is 

(2 (F + 47 r' 


X 


2 (F + 47t\ 

= V 3 C0S l^-j 


(33) 


(34) 


(35) 


where the quantity F is defined as 



F = f arccos ? < F < tt 


(36) 


It follows from equation (36) that 


or 




0.272166 


(37) 

(38) 
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and from equation (35) that 


X < 


[2 ( 57 t ' 


V 3 C0S V 


V3 


= 0.408248 


(39) 


These restrictions on the value of T and X should cause no difficulties for most applications. 


The value of S is obtained by eliminating the expression (EX 3 — R{) between equations (29) 
and (30), giving 

- X 


3X 2 V 1 — X 

and from equation (29), R\ is found to be given by 


(40) 


Rl = EX 3 + V 1 - X 




(41) 


The blending curve is defined uniquely by the quantity T Given T, the value of X is given by 
equations (35) and (36) and the values of E and R[ are given by equations (40) and (41), respectively. 

If the aircraft is in the flightpath segment defined by the blending curve, then the distance s(x p ) 
along the blending curve from the point 0 (fig. 5) to the foot of the perpendicular (x p , y p ) and the 
perpendicular distance 6y (fig. 5) must be determined. These distances are needed to determine the 
range and lateral deviation of the aircraft for guidance purposes. 

The equation of a line from the aircraft ( x a ,ya ) perpendicular to the blending curve, y = Ex 3 , is 
given by 

y-Va = 7*72 - Xa) (42) 

and, therefore, the coordinate x p at the intersection of this perpendicular and the blending curve is given 
by 




{x p - x a ) 


or 


3E 2 Xp - 3 Ey a x p + x p - x a = 0 


(43) 


An approximate solution of equation (43) may be obtained by setting x p = x a 
expression for x p in equation (43), and neglecting powers of i greater than one. 
the result is 

l = 3Exf(y a - Ex 3 ) 

1 - 3Sx a (2y a - 5Sx 3 ) 


4- 1, substituting this 
After simplification, 

(44) 



The coordinate y p is given by 


( 45 ) 



and the normalized perpendicular distance from the aircraft to the blending curve, or lateral deviation, 
6y, is given by 

<5 y = \J (x P - y a ) 2 + Wp - y a ) 2 (46) 

The normalized distance s(x) along the blending curve is obtained from the differential equation 

"2 


= ./i . /#V 
dx V \dxj 


(47) 


where y = E£ 3 . The formal solution of equation (47) is 

Xp 

s(x p ) = /(l 4- 9E*x^/ 2 dx 


(48) 


This integral cannot be evaluated in terms of the standard elementary functions; an approximate evalu- 
ation can be obtained by expressing the integrand in the power series 


(1 + 9S 2 ii 4 ) 1 / 2 = 1 + ^E 2 # 4 - ^E 4 x 8 

2 8 


which converges provided 


9S 2 S 4 < 1 


(49) 


(50) 


The largest value of x along the blending curve is X. Therefore, the limiting value of X for which the 
power series expansion converges is given by (eq. 40)) 


or 


-=■2 

9E 2 X 4 = —^2 < 1 
1-X 


X<-j= = 0.707107 

V2 


(51) 


(52) 


Since the solution of equation (33) is valid only for X < 0.408248, it follows that the power series 
expansion converges for the solutions of interest. 


Substituting equation (49) into equation (48) and performing the integrations gives 

, \ a | 9 ^2 5 9 A4 . 9 , 

s(x p ) = x p + —- 4 - Xp + 


(53) 


The maximum value of s(x p ) occurs when x p — X. Using the first three terms of the series gives 

s(X) * X(1 + ~E 2 X 4 - jE 4 * 8 ) (54) 
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and substituting for 2 from equation (40) gives 


s(X) « X 



( 55 ) 


Since the largest value of X is 0.408248, and the corresponding value of X is 0.680414, it follows 
from equation (55) that the the largest value of s(x p ) is 0.693643. The corresponding value obtained 
from equation (48) (with x p = 0.408248) by numerical integration is 0.693668. Therefore, the error in 
using the first three terms of equation (53) is always less than 0.000025. For a nominal turn radius of 
4,000 ft, the distance error is always less than 0.1 ft. 


ACQUIRING CURVE 


The purpose of the “acquiring curve” is to provide a smooth capture of the initial straight segment 
of the synthesized flightpath so that the pilot is relieved of the task of precisely aligning the heading 
of the aircraft’s velocity vector with that of the initial straight segment. By appending such a curve, 
the flightpath symbol on the HUD is made coincident with the ghost-aircraft symbol immediately at 
guidance select, and lateral guidance along the acquiring curve is provided through the ghost symbol to 
remove any initial misalignment between the velocity vector and the initial, straight flightpath segment. 

The first task in providing an acquiring curve and in determining those properties of it needed 
to drive the ghost symbol is to choose an appropriate analytical form. This form must be capable of 
representing a flightpath similar to one that the pilot himself would use in correcting the heading of 
the velocity vector; it must also be analytically tractable. Of the many reasonable possibilities, the one 
selected herein is a two-segment sine function to represent the second derivative of the acquiring curve 
(fig. 6). The rationale for this choice is based on the observation that, in correcting for a velocity- 
vector misalignment, the pilot would need to roll the aircraft toward the required flightpath track to 
produce a lateral acceleration relative to the track and then to roll in the reverse direction to reverse 
the lateral acceleration until the lateral velocity relative to the track is zero (“S” turn maneuver). It 
follows that, since the path curvature is proportional to lateral acceleration and since path curvature is 
approximately equal to the second derivative of the path, the form shown in figure 6 is a reasonable 
representation of what the second derivative of the path would look like if the pilot himself were to 
perform the misaligment correction task. A second important consideration in the selection of the form 
of the acquiring curve is that the curvature should be continuous as required to avoid discontinuities in 
the angle between the directions of flight of the ghost aircraft and the real aircraft. 

The equations for the form of the second derivative of the acquiring curve along with the first and 
second integrals, giving the slope and ordinate, are 


Segment 1 


0 < x < L\ 



A (l) L l f X \ 

= _-U_ C0S ^_J 


+ B (i) 


(56) 

(57) 
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Segment 2 


where 

L 
L i 

A (1)’ A (2) 



(?/)l = 


l (l) 


^1 • 
SH17T 


7T Z 



+ B^x + 


L\ < x < L 


(58) 



length of the acquiring curve measured along the x-axis (fig. 6) 
distance along the acquiring curve, measured along the x-axis, where the angle of 
bank changes sign (fig. 6) 
constants 
constants 


constants 


18 


4 



The six constants A( X ), A( 2 ), <7(1), and <7( 2 ) are set by the boundary conditions 


=E q -Vi 

\dxj i(x=0) 

(62) 

(¥) =0 

\dx J 2(x—L) 

(63) 

(dy\ _ 

Vdar/i(x=Li) \dx/ 2 (x=Li) 

(64) 

(?/)l(x=0) = 0 

(65) 

(^)l(x=Xi) = ( y)2[x=L x ) 

(66) 

(y) 2 {x=L) - 0 

(67) 


where 

E q heading of the HUD flightpath symbol at flightpath select 

heading of the initial flightpath segment 


Substituting the slopes and ordinates from equations (57), (58), (60), and (61) into equations (62)- 
(67) gives the following linear matrix equation whose solution vector contains the integration constants 


0 1 

ikzM 1 

(£-%) 0 

7 r 

0 Li 
0 0 
0 0 

The solution of equation (68) is 


0 0 0 ^ 


( A (U \ 


E q - 

-10 0 


A {2) 


0 

1 0 0 


%) 


0 

-£l 0 1 


B (2) 


0 

L 0 1 


C (l) 


0 

0 10; 


C (2) J 


0 / 


/ (Eg-WjML+Li) \ 


_£l 

7T 

k 1 

7 r 

0 

0 

0 

0 


A (i) 

A (2) 

B {i) 

B (2) 

C (l) 

C (2) 


\ 


J 


2L{L—L\) 

(Ea-VML-Lx) 
2L'~ 

(Ea-VJL 1 
2 1 

0 




(68) 


(69) 


Substituting the values of ^(i) , A (2) > -^(1) » -^(2) > Cfi), and C( 2 ) from equation (69) into equations (56)- 
(61) and expressing the result in terms of the nondimensional variables 
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~ x * _ y ? _ L i 

X L V L 1 L 

gives the final equations for the acquiring curve and its first and second derivatives: 


Segment 1 0 < x < L\ 

(<Py\ 

\dx 2 J i 


(Eg - gjMl + Lj) 


■m 


( i) i= (^)[ (1+£l)cos ,(i) +(1 _ £l) ] 

(£)l = sinw(£-) 4 (1 - h)x 


Segment 2 L\ < x < 1 


d 2 y\ __ (E q -^i)irLi (x-L\\ 

dx 2 ) 2 ~ 2(1 - Li) VI -Li) 

($) = J-^^ \L 1 cosJ X -^) + L 1 } 

\dx) 2 2 L \1- L\) J 

m = + £ l4 - u 


If (E q - is small, the curvature of the acquiring curve is approximately its second derivative 
and this curvature is proportional to the bank angle. The value of L\ should be such that the maximum 
and minimum bank angles are numerically equal in order to minimize the maximum bank angle needed 
during the acquiring segment of the flightpath. This result can be achieved by setting 


(d 2 y\ ^ = (<Py\ | 

\dx 2 Jix=h- \dx 2 ) 2\x=(^ 


Wdx 2 J i\x=Q I \dx 2 J 2 \x={^£±) 

Substituting for the second derivatives from equations (70) and (73) into equation (76) gives 

(1 + Li)tt _ Liir 

T x " 1-Li 


The solution of equation (77) is 


h = -= = 0.70711 
V2 


With the value of L\ given by equation (78), the maximum bank angle during the acquiring maneuver 
is given by 


<Pmax — arctan 


0.002055(£ 9 - %)vf 


where E q and W* are in degrees. For example, if vt = 200 ft/sec, (E q — = 10 deg, and L = 4,000 ft, 

then equation (79) gives <t> m ax = 11.6 deg. 
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CONCLUDING REMARKS 


Several problems with the HUD drive laws given in references 1 and 2 have been identified and 
solutions have been proposed. All these solutions have been flight tested successfully in the VSRA 
Phase IIA flight test program. Of all the problems addressed, the one whose solution appears to 
be marginal is that of acquiring the first segment of the synthesized flightpath when the aircraft’s 
course differs significantly from the heading of the first straight segment (identified as problem (7) in 
the Introduction). The synthesized-acquiring-curve approach given herein can be regarded only as an 
interim solution satisfactory for relatively small course errors (less than 15 deg). One pilot considered 
the “S” turn maneuver associated with the acquiring curve, coupled with the sensitivity of the ghost 
aircraft to lateral flightpath errors, to be impractical as an instrument procedure. A much better approach, 
although one requiring a more complicated algorithm, is to implement a two-circle flightpath synthesis 
requiring the pilot to correct his course to that of the first straight segment by making a simple turn in 
one direction only. 
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APPENDIX A 

FLIGHTPATH SYNTHESIS 


Preliminaries 

It is assumed that the onboard navigation sensors provide sufficient information and that there is 
sufficient processing capability to permit onboard continuous synthesis of a reference flightpath that is 
a member of a predefined class of flightpaths. When the synthesized flightpath is acceptable to the pilot 
he/she selects it by “freezing” the synthesizing process. 

The class of approach flightpaths used in conjunction with the HUD is illustrated in figure A-l. 
A set of constant parameters, identified in the following section, must be provided to the synthesizing 
system so that it can construct an appropriate flightpath within the predefined class of flightpaths. This 
flightpath is constructed as follows: from any aircraft position, the nominal, instantaneous, reference 
flightpath is defined as the line tangential to one of the circles shown in figure A-l and passing through 
the aircraft, followed by a segment of the circle and finally by a predefined straight-line segment. The 
flightpath ends either at touchdown or at the station-keeping point. The horizontal plane is divided into 
two parts by a line coincident with the final straight segment. This division determines which of the 
two circles will be used. If the aircraft is flown across the boundary during the process of continuous 
flightpath synthesis, the first segment of the flightpath moves from being tangential to one circle to 



Figure A-l. Class of reference flightpaths. 
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being tangential to the other. If the aircraft passes within either of the circles before the pilot selects a 
flightpath, the first segment switches from the circle entered to the other circle (fig. A-l). 


Required Data 


The following information must be provided to the system: 


S x , S y touchdown, or station-keeping point coordinates in the landing-pad reference frame 
landing-pad heading (0 < < 360 deg) 

f heading of the final flightpath segment (0 < <360 deg) 

D c nominal length of the final flightpath segment (fig. A-l) 

Rl( m in) nominal minimum radius of the horizontal flightpath circles 

Tp normalized nominal maximum length of the portion of the ordinate of the blending curve 

replacing a portion of a straight segment 
N 0 (min) minimum permitted value of uq 

ADac minimum desirable distance along the curved segment of the flightpath 


These parameters are generally constant and they have suitable default values, but they may be changed 
by the pilot any time prior to the time the reference flightpath is selected. After the reference flightpath 
has been selected, no further changes in these parameters can be made. In addition to these parameters, 
the navigation system must continuously provide the coordinates {x a , y a ) of the aircraft in the landing- 
pad reference frame. 


Synthesis 

Calculations before flightpath select- Figure A-2 shows the geometry of the flightpath. The 
coordinates of point C and the centers of the circles can be calculated from the basic data. The 
coordinates of point C are given by the equations 


Cx — Sx ~ Dc cos j 

(A-l) 

Cy = Sy — Dc sin A\&y 

(A-2) 

where A\D f is the heading of the final flightpath segment relative to the 
is given by the equation 

A = 0 < AW f < 360 deg 

x-axis of the landing pad. It 

(A-3) 

The coordinates S^ min ^ y ) of the centers of the minimum-radius circles are given by 

the equations 

R \ (L)x = C * - ' R l(min) 

(A-4) 

R< ulmn)y = Cy + l Rl(min) cos A j 

(A-5) 
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X 



where (z = +1) or (z = -1) denotes the circle for flight in the clockwise or counterclockwise directions, 
respectively. 

The appropriate circle depends on the distances lc \ i = +1, or z = —1, of the aircraft from the 
centers of the circles, given by 

W = J(xa- fl|;Ln)x) 2 + (*• - b 1(L in)/ where i = +1,-1 (A-6) 

with the circle-selection conditions 

if li +1) < Ri(rnin), then z = -1 
if z£ -1) < Ri( min y then % = +1 

if 4 +1) > Ri(min) 311(1 l c~ 1} > Rl(min ) and z c +1) > A~ l \ then i = -1 

if i +1) > Rl{min) 4 _1) > Rl(min ) 311(1 4 +1) < 4 _1) , then z = +1 

During the testing of the flightpath synthesis given in reference 1, the pilot was directed to make 
the heading change too abruptly when the heading change due to the turn from the initial to the final 

straight flightpath segments was less than about ±15 deg; the result was that large bank angles were 
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required for only short periods of time. This type of maneuver is an abnormal piloting technique even 
for a visual approach. The method used herein to overcome this problem is to allow the flightpath-circle 
radius to increase when the required heading change is less than about ±15 deg. This radius increase 
may be specified in many ways; the one adopted herein was chosen for its computational simplicity. It 
requires the calculation of the length l ac of the line AC and the angle 9q between the line AC and the 
final straight segment (fig. A-2). The value of l ac is given by 


lac — \J ( x a ~ Cx)^ + {ya ~ Cy)^ (A-7) 

The angle Dq is determined by first calculating the heading of the line AC, given by the equation 

-ipac = arctan ( ^ + W 0 < ip ac < 360 deg (A-8) 

VO# — x a / 

and then subtracting the heading of the final straight segment, taking into account the particular flightpath 
circle (value of i) being used. The calculation may be performed by using the equation 


i/q = 180(i + 1) - i(ipac ~ ^ f) 0 < vq < 360 deg 


(A-9) 


The value of uq must be limited in order to avoid excessive circle radii. This limited vq, designated 
%, is given by 


f ^0 



N 0{min) 


if J'o > N 0{min) 
if Uq < Ao(min) 


(A- 10) 


The radius of the flightpath circle is increased from the minimum, Ri( m i n y when the headings 
of the initial and final flightpath segments differ by less than about ±15 deg to a maximum, R ^ max ), 
whose value is controlled by specifying the minimum acceptable distance, A D ac , to be flown along the 
curved flightpath segment. The instantaneous circle radius, rj, and the blend parameter, v, are defined 
by the following relationships, which are shown graphically in figure A-3. 


If Uqi < 180 ADac/^Ri (min ) , then 

180AD ac 

n = — — 

™ 01 


(A- 11) 


The value of r/ obtained from equation (A- 11) must be bounded above to avoid, for as long as possible, 
flying into the interior of the circle prior to flightpath select— an event that would cause the synthesis 
procedure to switch to the other circle — and bounded below by the minimum radius. Specifically, 


r f 180 (lac - Rl(min)%)/™01 if r l > ^0(l ac ~ Ri {min) f p ) / nu oi 
1 l(min ) ^ r l ^ R l{min ) 

and the corresponding value of v is given by 

v = %Ri(min)/ r ll (A- 13) 
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Rightpath 

circle 



From inspection of the figure, the condition for the aircraft to be either outside or at the start of 
the blending curve is 

^ abc — 180r^o/7T + Ri(jnin)^p 
but for small values of vq., l a & c « l a c and uq^uq. Therefore, 

n < 180 (l ac - i?I(mm) T p)/^0 

Figure A-3 Geometrical condition restricting the flightpath radius. 


On the other hand, if uqi > 180AD ac / then 


r ll = R l(min) (A- 14) 

V = Tp (A-15) 


It follows from equations (A- 10) and (A-l 1) that the maximum possible flightpath circle radius, ify maa; ). 


is given by 


R 


180 AZ?ac 


(A- 16) 


Having determined the appropriate circle (value of i ) along with its radius, r^, and the blend 
parameter, v, the coordinates (r// x , ru y ) of this circle in the landing-pad reference frame are calculated 


i 
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by using the equations 


r llx = C X - ir H sin A^/ 
r lly = Cy + iru cos AVf 


(A- 17) 
(A- 18) 


Further calculations require knowledge of the angles 8ipi and A (fig. A-2). Angle Sipi is the angle 
between the line from the aircraft to the center of the appropriate circle (AL or AR of fig. A-2) and the 
x-axis of the landing pad. Angle A is the angle between AL or AR and the first flightpath segment. 
These two quantities are calculated from the equations 

8ipi == arctan ( — — ) 0 < 6tpi < 360 deg (A- 19) 

n \ r llx ~ x aJ 

A = — arcsin f^) (A-20) 

7 r \l c J 

where l c is the distance from the aircraft to the center of the circle that has been identified by the 
previous calculations; it is given by the equation 

lc = yj (xa ~ nix? + (Va ~ r U y) 2 (A-21) 


The track angle of the initial flightpath segment relative to the x-axis of the landing pad, may 
now be calculated from the equation 


6ipi — S'tpi — iX 0 < Sipi < 360 deg 
and the initial flightpath heading is given by 


^ + Stpi 0 <ipi< 360 deg 


(A-22) 


(A-23) 


The coordinates b x , b y of point B at the start of the circular segment may now be calculated using 
the equations 


bx = nix + lT ll sin (A-24) 

b v = r lly ~ ir ll cos b i>i (A-25) 


The angle vq that the circular segment subtends at the center of the circle is obtained using the 
equation 

uq = 180(i + 1) — — tyf) 0 < uq < 360 deg (A-26) 

The quantities i, b x ,b y , nh r llx^ r u y , and u 0 are calculated continuously up to the instant 

that the pilot selects the reference flightpath. After flightpath selection, these quantities remain constant 
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at their values at the instant of flightpath selection and they are subsequently designated by the uppercase 
symbols I, A W it B x ,B y , Ru, Ru x , R Uy , and N 0 . 


Calculations after flightpath select- The blending-curve constants are calculated immediately 
after flightpath selection. It should be recognized that there is a maximum possible value for v that 
may be less than the value of T p requested. This maximum is defined by two possibilities: first, the 
appropriate solution resulting in a proper fit of the cubic blending curve to both the straight and circular 
segments holds only if 


v < 



(A-27) 


and second, the value of v cannot be greater than that which causes the blending curves for entry and 
exit to meet in the middle of the circular segment. The condition expressing this limitation is 


v < sin 


No 

2 



N 0 


< 2 arcsin 



(A-28) 


The final value of v is then given by the lesser of the limiting values given by equations (A- 13) 
or (A-15), (A-27) and (A-28). This value is designated T 

The blending-curve constants are obtained by first defining a quantity F as follows: 



arccos 



(A-29) 


The normalized distance X measured along the straight segment over which the blending curve 
replaces the circular arc segment is then given by 


and the normalized distance X measured along the straight segment is then given by 

x = x + r 


(A-30) 


(A-31) 


The defining constant of the cubic blending curve, E, is then calculated from the equation 


X 

3X2\/l-X 2 


(A-32) 


If T = 0, then X and X are both zero and S from equation (A-32) is indeterminate. In order to recover 
the flightpath-synthesis equations of reference 1, which are appropriate to the case where T = 0, it is 
necessary to define E to be zero. 
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The true radius Ru of the flightpath circle is obtained from the nominal radius Ru by using the equation 

Ril R “, (A-33) 

zx 3 + yi-x 

If the quantity E is defined as follows, 

E = f arcsin X (A-34) 

then the true flightpath circular segment subtends an angle Nq given by 

N 0 = N 0 - 2 E (A-35) 

Relationships (A-29)-(A-35) complete the calculation of the blending-curve constants; we now will 
calculate the key ranges (measured along the true reference flightpath) at the points C\, B\, B 2 , and 

A. First, it is necessary to calculate the lengths of the two principal segments of the blending curve 
corresponding to the normalized ordinates T and X These lengths can be obtained from the equation 



s(x) =x + 0.9H 2 # - 1, 125H 4 x 9 

(A-36) 

where 



X 

s(x) 

normalized blending-curve abscissa 
normalized length of the blending curve 


The key 

ranges D cl ,D C2 ,D bl , and D b<2 are then given by 



D c 1 = D c - Ruf 

(A-37) 


D C2 = D C1 + Rus(X) 

(A-38) 


n n , vRllNo 

0^=002+ lg0 

(A-39) 


D bl = D b2 + D c 2 — D Cl 

(A-40) 

The instantaneous range D a of the aircraft at flightpath selection is given by 



D a = l c cos A - RuT + D bl 

(A-41) 


The quantities R U ,E, N 0 ,T , D Cl , D C2 , D bl , D b2 , and D a are calculated once only for subsequent use 
in the guidance equations (appendix B). 
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APPENDIX B 
GUIDANCE 


Lateral Guidance 


Lateral guidance is provided by the ghost-aircraft symbol, and the lateral-guidance law is provided 
by the geometric relationship between the flightpath symbol and the ghost (fig. 16 of ref. 1). Calculation 
of the azimuth angle of the ghost, n g , requires knowledge of the range, d, the lateral offset of the aircraft 
from the reference-flightpath track, 6y, and the heading of the reference flightpath, ipt. In addition, the 
ghost roll angle, 4> g , given by 


180 

4> g = — arctan 



(B-l) 


requires knowledge of the reference-flightpath curvature, 1/p, at the aircraft’s position (identified by its 
range, d). 


The quantities I , Ad/^, Wj, Bx, By, Rn, Riix> Rlly> Aqj 1, -Dei* Dc <%, D b ^, nnd Da have 
been calculated by the flightpath-synthesis algorithm (appendix A). The reference flightpath is now 
modified slightly by the addition of the acquiring segment; this segment is defined by the following five 
quantities: 

L overall length of the acquiring curve segment 

Kdb ratio of the lengths of the acquiring curve and the initial straight segment 
L\ normalized length from the start of the acquiring curve to the point of roll reversal 

E q heading of the HUD flightpath symbol at the instant of flightpath select 

heading of the initial flightpath segment 


The quantity L\ is normally set equal to l/\/2, but it may be changed by the pilot. The length L 
is normally set equal to the distance of the aircraft from the start of the turn entry, D a — D bl , but it 
may be changed to a fraction of this distance, K ab , by the pilot. Specifically, 


L = K ab (D a -D h ) (B-2) 

The heading E q (0 < E q <360 deg) is set equal to e q 4- \& a (e q is given by eq. (15) of ref. 1) at the 
instant of flightpath select and is constant thereafter. The equations needed to determine the quantities d, 
6y, and ipt depend on which segment of the reference flightpath corresponds with the aircraft’s location 
(fig. B-l). 


For segments AA\, A\B\, B 1 B 2 , and B 2 C 2 , the normalized coordinates of the aircraft (X 2 a , y2a) 
in an axis system with origin at B and a>axis along the line AB are required and are given by 


%2a = [{x a ~ B x ) cos A % + (y a - B y ) sin A '& l ]/Rn + f (B-3) 

V2 a = -I[{xa ~ Bx) sin A^ i - ( y a - B y ) cos A ^i]/Ru (B-4) 
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Also required for calculations in segments AA\ and A\B is the range, D ai , of the point A\ at the end 
of the acquiring curve, given by 

D ai =D a -L (B-5) 

The parameters within the acquiring segment AA\ defined by D a > d > D ai 

where 


are given by 


if 0 < x < L\ 


d = D bl - Riix 2a 
x = (D a — d)/L 


6y - IRuma 
7rW 


7 TL^ 


VI 


360 


Li(l + Li) . (x 
- — — - - sm7r — 

\L X 


ix 


+ (1 - Li)x 


if>t = + 


vi 


360 


, A . . 


1 

P 


7T 2 (l+£l)fr 

360LLi 


sin 


y A . 

■© 


(B-6) 

(B-7) 

(B-8) 

(B-9) 

(B-10) 
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and by 


where 


Sy IR ^ 2a 360 L 

TT^ 


nL^, f Xl(l-£l) sjn , (x-U'j + L 1 i-L 1 




VI 


360 

1 _ -it 2 Li'S. 

P ~ 360L(1 - Li) 


9 /x-L 1 \ 9 ' 

Licos7r( — + L\ 

\1 — Li / J 


v% 

— — sm 


. fx-Li\ 
m 7r ( — 

\l-UJ 


(B-ll) 

(B-12) 

OB-13). 


%i = E q -Vi (0 < < 360 deg) 


(B-14) 


if L\ < x < 1. 


It should be noted here that if an acquiring curve is not required {K ab = 0), then the calculations of 
equations (B-5)-(B-14) must be omitted. 

The parameters within the initial straight segment A\B\ defined by D ai > d > D bl 


where 


d = Dbi ~ Rll%2a 

(B-15) 

are 


6y = IRnha 

(B-16) 


#- 

II 

•e 

«s>. 

CB-17) 


1=0 

(B-18) 


P 


The parameters within the turn-entry blend segment B\B 2 defined by D bl >d> D b2 
where 

d =D bl - Riis(x p ) 


(B-19) 


are 


8y 

1 

P 


[sign ( y 2a ~ S xla)lI R liy/i^2a ~ *p) 2 + (V2a ~ tip? 

~ —— arctan (35 x 2 ) 

7 r ■■ 

I6Ex p 

Wi I + 9=2iJ)3/2 


(B-20) 

(B-21) 

(B-22) 
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For the segment B 2 C 2 , the coordinates of the aircraft (xi a ,yi a ) are required in an axis system 
whose origin is at the center of the appropriate circle and whose z-axis is parallel to the final flightpath 
segment (fig. A-2). These coordinates are given by 



xu = (x a - Riix ) cos A % + (y a - Ru y ) sin AW f 

(B-23) 


Via = ~ Rux ) sin A W/ - ( y a - R Uy ) cos AW f ) 

(B-24) 

In addition, the angle u subtended at the center of the circle by the remaining part of the circular-arc 
segment is required and is given by 


v — arctanf — — ^ 0 < v < 360 deg 

7T V Via / 

(B-25) 

The parameters within the circular segment B 2 C 2 defined by > d > D C2 


where 

d - D C2 + 18 Q Ru(v E) 

(B-26) 

are 

la = \J ( X a - Rllx ) 2 - ( Va ~ Rlly) 2 

(B-27) 


Sy = I{Ru - la) 

(B-28) 


1 

II 

(B-29) 


1 _ I 
p Ru 

(B-30) 


For segments C 2 C 1 and C\S, the normalized coordinates of the aircraft (£ 2 a » V2 a) are redefined 
relative to an axis system with origin at C\ and x-axis along the line SC; they are now given by 

£ 2 a = -[(*0 - S x ) cos A + (y a - S y ) sin AW/ + D Cl ]/Ru (B-31) 

V2 a = -I[-(ya ~ Sy) cos A W / 4- (x a - S x ) sin AW f]/Ru (B-32) 

The parameters within the tum-exit blend segment C 2 C 1 defined by D C2 >d> D Cl 


where 

are 


d = D Cl + Ru$(x p ) 


Sy 

i’t 

1 

P 


[sign(y2a - Zx2a)}IRll^/(x2a ~ x p ) 2 + {y 2a - Vp ) 2 

W t — arctan (3Sa£) 

J ft F 

IQExp 

Bui l + 9S2£4)3/2 


(B-33) 

(B-34) 

(B-35) 

(B-36) 
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The parameters within the final straight segment C\S defined by D Cl > d 


where 


= D c 1 + R U X 2a 

(B-37) 

6y = I Rama 

(B-38) 

$ 

II 

■e 

(B-39) 

1 = 0 

(B-40) 


P 


The quantities x p and y p are the normalized coordinates of the foot of the perpendicular from the 
aircraft to the blending curve; they are given by 


x p — X2a 


1 + 


3^^2a(^2a ,=, ^2 a ) 


1 — 3c.X2a(2y2a — ^"®2a) 


Vp 



(B-41) 

(B-42) 


The quantity s(x p ) is the normalized distance of the point (x p ,y p ), measured along the blending 
curve from the point where the transition curve and a straight segment meet; it is given by 

s(x p ) — x p + 0.9 E^Xp - 1.125 E^x p (B-43) 


Vertical Guidance 

Vertical guidance, like lateral guidance, is provided by the ghost-aircraft symbol. At the instant of 
guidance select, the vertical flightpath is defined as shown in figure B-2. This flightpath consists of a 
constant-altitude segment, AJ\, at the altitude of flightpath select, an approach segment, J 2 S, with a 
constant flightpath angle, T, and a circular connecting segment, J\ J 2 The following information must 
be provided to the system in order to define this flightpath: 

H r i aircraft altitude at guidance select 

H hover altitude 

Ry radius of the vertical flightpath circle, J\ 3 2 
T approach flightpath angle 

The quantities H, Ry, and T are set by the pilot, with default values of 50 ft, 30,000 ft, and -3 deg, 

respectively. 
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From figure B-2, the ranges at points J, J\, and J 2 are given by 


D j = (H-H ri )/um( r) 

(B-44) 

D h = D i ~ -Rotand ) 

(B-45) 

Dj 2 = Dfr -F Rv tan(r) 

(B-46) 


The elevation angle of the ghost relative to the aircraft, rjg, is expressed in reference 1 (eq. (19)) 


as 


Vg 


180 (hg — h\ 

= — arctan -f- 

7T V 6x gh ) 


(B-47) 


where 

h g altitude of the ghost aircraft 

h altitude of the aircraft 

8x g h lead distance of the ghost for vertical tracking 


The ghost altitude is given by 

hg = hr + hXgfokyjg tUO 

where 

h r aircraft reference altitude 

7 1 reference flightpath angle at the aircraft’s range 

k wg ghost lead blending gain 


The ghost lead blending gain k wg is defined in reference 1 (eq. (22)). 


(B-48) 
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From equations (B-47) and (B-48), r\g may be expressed as 


% 


— arctan 


(k wg tan7t - 


8h \ 
^•Egh ' 


where 8h = h - h r is the altitude error of the aircraft relative to the reference flightpath. 


(B-49) 


Determination of rj g , therefore, requires knowledge of the altitude error, 8h, and the reference 
flightpath angle, 7 f. These quantities depend on which segment of the reference flightpath corresponds 
with the aircraft’s location. 


The parameters within the constant-altitude segment AJ\ defined by D a > d> Dj x 
are 


It = 0 (B-50) 

8h — h- H ri (B-51) 

The parameters within the circular segment J\ J 2 defined by Dj 1 >d> Dj 2 

are 


(B-52) 
(B-53) 

The parameters within the constant-flightpath-angle segment J 2 S defined by Dj 2 > d > 0 
are 


It = arcsin[(d - Dj^/Rv] 

8h — h — H r i + Rv(l — cost 1 ) 


lt = T (B-54) 

8h = h-H + d tan(T) (B-55) 
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